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Abstract 

In this paper, we want to understand the Proudman resonance. It is a resonant 
respond in shallow waters of a water body on a traveling atmospheric disturbance 
when the speed of the disturbance is close to the typical water wave velocity. We 
show here that the same kind of resonance exists for landslide tsunamis and we 
propose a mathematical approach to investigate these phenomena based on the 
derivation, justification and analysis of relevant asymptotic models. This approach 
allows us to investigate more complex phenomena that are not dealt with in the 
physics literature such as the influence of a variable bottom or the generalization of 
the Proudman resonance in deeper waters. First, we prove a local well-posedness 
of the water waves equations with a moving bottom and a non constant pressure 
at the surface taking into account the dependence of small physical parameters and 
we show that these equations are a Hamiltonian system (which extend the result of 
Zakharov [3S1)- Then, we justify some linear asymptotic models in order to study the 
Proudman resonance and submarine landslide tsunamis; we study the linear water 
waves equations and dispersion estimates allow us to investigate the amplitude of the 
sea level. To complete these asymptotic models, we add some numerical simulations. 


1 Introduction 

1.1 Presentation of the problem 

A tsunami is popularly an elevation of the sea level due to an earthquake. However, 
tsunamis induced by seismic sources represent only 80 % of the tsunamis. 6% are due 
to landslides and 3% to meteorological effects (see the book of B. Levin and M. Nosov 
|2nj ). Big traveling storms for instance can give energy to the sea and lead to an elevation 
of the surface. In some cases, this amplification is important and this phenomenon is 
called the Proudman resonance in the physics literature. Similarly, submarine landslides 
can significantly increase the level of the sea and we talk about landslide tsunamis. In 
this paper, we study mathematically these two phenomena. We model the sea by an 
irrotational and incompressible ideal fluid bounded from below by the seabed and from 
above by a free surface. We suppose that the seabed and the surface are graphs above 
the still water level. We model an underwater landslide by a moving seabed (moving 

‘benjamin, melinand@math.u-bordeaux.fr 


1 



bottom) and the meteorological effects by a non constant pressure at the surface (air- 
pressure disturbance). Therefore, we suppose that b{t,X) = bo{X) + bm{t,X), where bo 
represents a fixed bottom and bm the variation of the bottom because of the landslide. 
Similarly, the pressure at the surface is of the form P+Pj-^f, where Pref is a constant which 
represents the pressure far from the meteorological disturbance, and P{t,X) models the 
meteorological disturbance (we assume that the pressure at the surface is known). We 
denote by d the horizontal dimension, which is equal to 1 or 2. X G stands for the 
horizontal variable and z G M is the vertical variable. H is the typical water depth. The 
water occupies a moving domain fit := {{X,z) G — H + b{t,X) < z < C(t,X)}. 

The water is homogeneous (constant density p), inviscid, irrotational with no surface 
tension. We denote by U the velocity and $ the velocity potential. We have U = Vx,z^- 
The law governing the irrotational fluids is the Bernoulli law 

dt^ + +gz=- (Pref - p) in Oi, (1) 

2 p 

where V is the pressure in the fluid domain. Changing <h if necessary, it is possible to 
assume that Pref = 0. Furthermore, the incompressibility of the fluid implies that 

Ax,z^ = 0 in (2) 

We suppose also that the fluid particles do not cross the bottom or the surface. We denote 
by n the unit normal vector, pointing upward and the upward normal derivative. 
Then, the boundary conditions are 

dtC - Vl + |VC|29n$ = 0 on {z = C{t, X)}, (3) 

and 

dtb - v^l + |V6|29n$ = 0 on {z = -H + bit, X)}. (4) 

In 1968, V. E. Zakharov (see [33]) showed that the water waves problem is a Hamiltonian 
system and that 'll), the trace of the velocity potential at the surface (ip = <I>| 2 =^), and 
the surface (p are canonical variables. Then, W. Craig, C. Sulem and P.L. Sulem (see 
ffH and m) formulate this remark into a system of two non local equations. We follow 
their construction to formulate our problem. Using the fact that satisfies Q and Q, 
we can characterize $ thanks to C and ip = 

f Ax ^<I> = 0 in Qt, 

\ ’ ^ _ (5) 

\ = p; , ,/lTWWdn%=-H+b = dtb. 

We decompose this equation in two parts, the surface contribution and the bottom 
contribution 


such that 
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( 6 ) 


Ax,= 0 in 

^ > V^+]Wdn^%=_H+b = 0 , 

and 

Ax,:<&^ - 0 in 

<!>';„( = 0 , Vl + |V6P8„<t'«,=-H+i, = Dll’- 

In the purpose of expressing Q with ( and ifj, we introduce two operators 
one is the Dirichlet-Neumann operator 

G[C, 6] : V’ ^ (8) 

where satisfies The second one is the Neumann-Neumann operator 

G^^iC, b] : dtb ^ Vl+W^dn^%=^, (9) 

where satisfies (0. Then, we can reformulate (0 as 

dtc-G[c,bm = G^^[c,b]idtb). ( 10 ) 

Furthermore thanks to the chain rule, we can express (Vx, 2 *^) 1 ^=^ and 

in terms of ifj, (, G[C, ^KV') G^^[C, b]{dtb). Then, we take the trace at the surface of 

Q (since there is no surface tension we have 'P\z=(; = P) and we obtain a system of two 
scalar equations that reduces to the standard Zakharov/Craig-Sulem formulation when 
dtb = 0 and P = 0, 


( 7 ) 

The first 


dtc-G[c,bm = G^^^^[c,b]{dtb), 

1 (G[C, 6 ]«-) + G™|C,I. 1 (W + VC-V^)" p 

aV + K+jIVV-l -2-jiXTvCp)-= ■/ 


( 11 ) 


In the following, we work with a nondimensionalized version of the water waves equations 


with small parameters £, (3 and /r (see section 2.1). The wellposedness of the water waves 
problem with a constant pressure and a fixed bottom was studied by many people. S. 
Wu proved it in the case of an infinite depth without nondimensionalization m and 
[32]). Then, D. Lannes treated the case of a finite bottom without nondimensionalization 
([T7]). T. Iguchi proved a local wellposedness result for ^ small enough in order to justify 
shallow water approximations for water waves ([IS]), and D. Lannes and B. Alvarez- 
Samaniego showed, in the case of the nondimensionalized equations, that we can find an 
existence time T = where Tq does not depend on e, (3 and /x (|6]). More recently, 

B. Mesognon-Gireau improved the result of D. Lannes and B. Alvarez-Samaniego and 
proved that if we add enough surface tension we can find an existence time T = — 
where Tq does not depend on e and /x ([23]). T. Iguchi studied the case of a moving 
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bottom in order to justify asymptotic models for tsunamis ([IS]). Finally, T. Alazard, 
N. Burq and C. Zuily study the optimal regularity for the initial data m) and more 
recently, T. Alazard, P. Baldi and D. Han-Kwan show that a well-chosen non constant 
external pressure can create any small amplitude two-dimensional gravity-capillary water 
waves (ID). We organize this paper in two part. Firstly in Section we prove two 
local existence theorems for the water waves problem with a moving bottom and a non 
constant pressure at the surface by differentiating and ” quasilinearizing” the water waves 
equations and we pay attention to the dependence of the time of existence and the size 
of the solution with respect to the parameters e, /?, A and /i. This theorem extends the 
result of T. Iguchi m ) and D. Lannes (Chapter 4 in (Hj). We also prove that the water 
waves problem can be viewed as a Hamiltonian system. Secondly in Section]^ we justify 
some linear asymptotic models and study the Proudman resonance. First, in Section [3T] 
we study the case of small topography variations in shallow waters, approximation used 
in the Physics literature to investigate the Proudman resonance; then in Section 3.2 we 


derive a model when the topography is not small in the shallow water approximation; 
and in Section 3.3 we study the linear water waves equations in order to extend the 
Proudman resonance in deep water with a small fixed topography. Finally, Appendix [A] 
contains results about the elliptic problem and Appendix contains results about 
the Dirichlet-Neumann and the Neumann-Neumann operators. Appendix comprises 
standard estimates that we use in this paper. 


1.2 Notations 

A good framework for the velocity in the Euler equations is the Sobolev spaces But 
we do not work with U but with the trace of and U = Vx, 2 ^. It will be too 
restrictive to take -0 in a Sobolev space. A good idea is to work with the Beppo Levi 
spaces (see m)- For s > 0, the Beppo Levi spaces are 

:= 1^/; E LL(M'^), VV' E . 

In this paper, C is a constant and for a function / in a normed space (A, I'l) or a 
parameter 7, C(|/|,7) is a constant depending on |/| and 7 whose exact value has non 
importance. The norm | • 1^2 is the L^-norm and | • |oo is the L°°-norm in Let 

/ E and m E N such that £ L°°(M'^). We define the Fourier multiplier 

f{D) : ^ L2(M'^) as 

Vn E fiD^m = fiOm- 

In we denote the gradient operator by V and in H or S = x (—1,0) the gradient 
operator is denoted Vj)c, 2 . Finally, we denote by A := y^l -|- |L)p with D = —iV. 
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2 Local existence of the water waves equations 


This part is devoted to the wellposedness of the water waves equations (Theorems 2.3 


and 2.4). We carefully study the dependence on the parameters e, /3, A and /i of the 
existence time and of the size of the solution. Contrary to m and we exhibit the 
nonlinearities of the water waves equations in order to obtain a better existence time. 


2.1 The model 

In this part, we present a nondimensionalized version of the water waves equations. In 
order to derive some asymptotic models to the water waves equations we introduce some 
dimensionless parameters linked to the physical scales of the system. The first one is the 
ratio between the typical free surface amplitude a and the water depth H. We define 
e := jj, called the nonlinearity parameter. The second one is the ratio between H and the 
characteristic horizontal scale L. We define fi called the shallowness parameter. 

The third one is the ratio between the order of bottom bathymetry amplitude Obott and 
H. We define /3 := called the bathymetric parameter. Finally, we denote by A the 
ratio of the typical landslide amplitude abott,m and Obott- We also nondimensionalize the 
variables and the unknowns. We introduce 


X' = — z' = — C = - b' = ^ h' = h' = t' = t 

■L/ II d Q^bott ^bott ^bott.m ^ 


= 


H 


aL^/gH 




- - tP' = — — p' = 

’ aL^/^ ’ apg' 


( 12 ) 


where 


0) = {{X', z') G - 1 + X') <z' < eC'it', X')}. 

Remark 2.1. It is worth noting that the nondimensionalization of , p and t comes 
from the linear wave theory (in shallow water regime, the characteristic speed is y/gH). 
See paragraph 1.3.2 in m- Let us explain the nondimensionalization of . Consider 
the linear case 


I = 0, - R < ^ < 0, 

\ ^%=o = 0 .= dtb. 

A straightforward computation gives If the typieal wavelength is 

L, the typical wave number is Furthermore, the typical order of magnitude of dtb is 
ahott,rnV9H _ ^ ^/jg of magnitude of in the shallow water case is 

^ y/9^^bott,m sinh(27r-^) \/ QHcl bott,mLI 

27r L cosh(27r^) L 
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For the sake of clarity, we omit the primes. We can now nondimensionalize the water 
waves problem. Using the notation 

^x,z ■= (\//^^x , Y and := /rA^ + 5^, 

the water waves equations 0 become in dimensionless form 


--^G^ [eC, / 56 ] (V^) = K, m {dth ), 

- (1 + .W) -^ 


In the following is the upward conormal derivative 

= ii-( ° 


-P. 

(13) 




0 ij 


Then, The Dirichlet-Neumann operator G^[eC,/36] is 


G^[eC, I3bm ■■= + = -frsVC • , 

where <I>'^ satisfies 

r = 0 in Ot , 

I ‘^iz=eC “ ^ ’ ^^^%=-l+/3b — 0) 
while the Neumann-Neumann operator (3b] is 


(14) 


(15) 


G^^[sC,(3bmb) := + + , (16) 

where <I>^ satisfies 


J = 0 in Sil . 

1 = 0 , s/lTMW9„>I>®|,._i+,, = dA 

Remark 2.2. We have nondimensionalized the Dirichlet-Neumann and the Neumann- 
Neumann operators as follows 


G[c,bm = 


aLy/gH 


G^[e( = 


^bott, 




G^^^[e(,^t;]{d,,lf). 


L 
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We add two classical assumptions. First, we assume some constraints on the nondi- 
mensionalized parameters and we suppose there exist pmax > 0 and ^max > 0, such 
that 


0 < £, /?, /3A ^ 1 , ^ ^ Pmax and // ^ /^max- (1^) 

Furthermore, we assume that the water depth is bounded from below by a positive 
constant 


3 ^min ^ 0 ? eC + l- Ph> 

^min* (19) 

In order to quasilinearize the water waves equations, we have to introduce the vertical 
speed at the surface w and horizontal speed at the surface (F. We define 


w := ^[eC) 



and 


G^K,/36](V^) +/i^G;y^K,/36](9i6) +£/iVC • 

1 + £V|VCP 


( 20 ) 


V := F[eC, lib] ^dtl^ =V^- ew[eC, pb] (^^1, ^dtb'^ VC- (21) 

2.2 Notations and statement of the main results 

In this paper, d = 1 or 2, to > ^, N ^ 'N and s > 0. The constant T > 0 represents 
a final time. The pressure P and the bottom b are given functions. We suppose that 
5 g and P G We denote by Mn a constant 

of the form 

Mn = C f - , /imax) £|CI Hmax(to+2,Ar) , /3|6|^oo^inax{to+2,iV) j . (^2) 

\ ^min t X J 

We denote by U := (C, V')* the unknowns of our problem. We want to express © as a 
quasilinear system. It is well-known that the good energy for the water waves problem 
is 


£^{u) = m\^^,+ Y. (iCHii 2 + iW(a)ii 2 ), (23) 

where C(a) := '4’{a) ■= “ £wd^( and ip := This energy is motivated 

by the linearization of the system around the rest state (see 4.1 in US]). ^ acts as the 
square root of the Dirichlet-Neumann operator (see [l9|). Here, C(o) and '^(q) are the 
Alinhac’s good unknowns of the system (see [I] and [3| in the case of the standard water 
waves problem). We define := (C(a)) V'(a))*- We can introduce an associated energy 
space. Considering a T > 0, 
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:= {U x , £^{U) E L“([0,r])}. (24) 


Our main results are the following theorems. We give two existence results. The first 
theorem extends the result of T.Iguchi (Theorem 2.4 in [16]) since we give a control of 
the dependence of the solution with respect to the parameters e, (3 and fj. and we add a 
non constant pressure at the surface and also extends the result of D.Lannes (Theorem 
4.16 in m), since we improve the regularity of the initial data and add a non constant 
pressure pressure at the surface and a moving bottom. Notice that we explain later 
what is Condition ( [^ (it corresponds to the positivity of the so called Rayleigh-Taylor 
coefficient). 

Theorem 2.3. Let A >0,to> N > max(l,to)+3, , b £ W^’°^{R+; {R‘^)) 

and P £ W^’°°{R+]H^+^{R^)) such that 


gN ^ ^ < A. 

are satisfied 
with initial 

data U^. Moreover, we have 


We suppose that the parameters e, fi, p, X satisfy ( |18[ ) and that (19) and (29) 
initially. Then, there exists T > 0 and a unique solution U £ ES to (13 


T = min 


To __ 

max(e,/3)’ ^ \dtb\i 


To 


+ I^Tlj 


0 r/iV 
iix , 


— = and sup £^{U) = c^, 
To ie[o,T] 


with cf — C ^^4, , ^max, Pmax, Hfi ’ \ Hfi ) ' 

Notice that if dfi) and P are of size max(e,/3), we find the same existence time that 
in Theorem 4.16 in m- The second result shows that it is possible to go beyond the 
time scale of the previous theorem; although the norm of the solution is not uniformly 
bounded in terms of e and fi, we are able to make this dependence precise. This theorem 
will be used to justify some of the asymptotic models derived in Section over large 
time scales when the pressure at the surface and the moving bottom are not supposed 
small. We introduce 6 := max(e,/3^). 

Theorem 2.4. Under the assumptions of the previous theorem, there exists To > 0 such 
that U £ E^^ . Moreover, for all a £ [O, , we have 

7J 


4 = c‘,sup £"((7) < c’ = c(a, f- 


|VP| 








Notice that when dtb and P are of size max(e, /?), the existence time of Theorem 2.3 


IS 


better than the one of Theorem 2.4 Theorem |2.4| is only useful when dtb and P are not 


small. Notice finally, that Condition (29) is satisfied if s is small enough. Hence, since 


in the following, e is small, it is reasonable to assume it. 
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2.3 Quasilinearization 

Firstly, we give some controls of and |^V’(q!)|h= with respect to the energy S^{U). 

Proposition 2.5. Let T > 0, to > ^ and N > 2 + max(l,to)- Consider U G , 
b G such that ( and b satisfy Condition (19) for all 0 < t <T. We 

assume also that yi satisfies (18). Then, for 0 < t < T, for a G N'^ with |a| < N — 1 and 
for 0 < s < N — 


+ mia)\m + m\H^ < Mn£^{UY 2 + illlMN\dth\L^H^ . 


/3A 


Proof. For the first inequality, we have thanks to Proposition C.l[ 


< IW(a)lL2 +e|q}(:«;a“C)lL2, 


< IW(a)ll,2 + — |uz5“C| 

/i4 


1 . 


But Tp G Then by Proposition B.8 w G 

Proposition |C.2t we obtain 


and a“C e 


Using 


< m^^)\L2+Ce 


w 




ICI//JV < |ipV'(a)lL2+MAre|C|iyiV 

m 


The other inequalities follow with the same arguments, see for instance Lemma 4.6 in 

dn. □ 

The following statement is a first step to the quasilinearization of the water waves equa¬ 
tions. It is essentially Proposition 4.5 in m and Lemma 6.2 in [16]. However, we 
improve the minimal regularity of U (we decrease the minimal value of to 4 in dimen¬ 
sion 1) and we provide the dependence in dtb which does not given in [T6|. For those 
reasons, we give a proof of this Proposition. 

Proposition 2.6. Let to > ^, T > 0, N > max(to, 1) -|- 3, 6 G and 

U G , such that ( and b satisfy Condition (19) for all 0 < t <T. We assume also 
that fj, satisfies (18). Then, for all a G N'^, 1 < |a| < N, we have. 


1 


A/3 


1 


/3A 


5“ ( -G^[eC, mW + —G^ =-GuH, m (V’(a)) + —C [<, f3b] {d^dtb) 

JX £ J T ^ 


^l-{|a|=7V}^ ■ (C(o)lA) -|- Ra‘ 


Furthermore R„ is controlled 


/3A 


\Ra\L^ < MN\i£C,mH^£ (U)-^ + —Mn\ dtb\L^H^ . 

£ t X 
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Proof. We adapt and follow the proof of Proposition 4.5 in [T9]. See also Proposition 


6.4 in |16) . Using Proposition B.13 we obtain 


1 


A/3 


1 


/3A 


£ / /^ ^ 


- £1{|„|=W}V • (C(«)Z) + (v • [d%v 

+ T?a, 


where U = Vfef, /35]('0, B) is defined in Equation ( |7^ and Ra is a sum of terms of the 
form (we adopt the notation of Remark B.12 in Appendix |B.3|) 


where j is an integer and ...,P and u are multi-index, and 

\i'\ + W\=N, 

i<i<j 

with (j, |z^|) / (IjA", 0) and (0,0, A"). Here is such that max|t*| = In 

i<*<i 

particular, 1 < < N. We distinguish several cases. 

a) + |z/| < N — 2 and |t*°| < N — 3 or + |i/| < N, < N — 3 and \i>\ < N — 2 : 


Applying the second point of Theorem 3.28 in m and the hrst point of Proposition 
with s = 2 and to = niin(to, |)) we get that 


B.15 




\^d-^\^, + ^\d'^dtb\L2 


and the result follows by Proposition 2.5 
b) |t^o| = _/v — 2 and |z^| = 0 ,1 or 2 : 


We apply the fourth point of Theorem 3.28 in m and the second point of Proposition 
B.15 with s = 2 and to = max(to, 1), 


\Aj,,,,\L.<MM{ed‘'\,pd^^b)\ sillied^ C,/ 39 ^ b)\H^ 


l^lo 


m2d''f;\HN-2 + ^\d''dtb\HN-2 


Then, we get the result thanks to Proposition 2.5 
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c) = i with |i| = — 1, |i^| = j = 1 : 


We proceed as in Proposition 4.5 in m, using Theorem 3.15 in m and Propositions 

Bi^iBTnra 


d) = N — I and |z^| = 0 : 


Here j = 2 and \r\ = 1. For instance we consider that = 1 and \r\ = 1. Using the 


second inequality of Proposition B.15 we have 


(dtb).C, C; 9^' b, || C, |^J ft) |^. 


Furthermore, using two times Proposition B.13 we get 




= -^dG^[eC,/36] .(a^'c,0) 

- ^G^K,/36] 

-eV- (9"'Cf^U(V^,0).(5""c,0)) 

+ (3dG^^[£C,(3b] (9"'6U(V’,0)) .(0,5^'5) 

+ PG^^^H,l5h] (a‘^6d£(V>,0).(0,9^'6)) . 


The control follows from the first inequality of Theorem 3.15 and Proposition 4.4 in m, 


and Propositions B.14 B.8 and B.ll 


e) \iy\ = N — 1 and |i^“| =1 : 


Here, j = 1- It is clear that 


—dG^^{d'^dtb).{d‘-\-,d‘-"b) 


L2 


< —Mat \ dtb\jjN 


Furthermore, 


-dG^{d'^^).id^\-,d^"b) = -dG^i^l^^,)).{d‘-\-,d‘-"b) + -^dGj^wd^c) ■{d‘'\-,d‘-"b). 

A* y \y / 

Then, using Theorem 3.15 in m and Proposition |2.5[ we get the result. □ 
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This Proposition enables to quasilinearize the first equation of the water waves equations. 
For the second equation, it is the purpose of the following proposition. 


Proposition 2.7. LetT > 0, > max(to,l) + 3, b G ■, andU G , 


such that C and b satisfy (19) for all 0 < t <T. 
Then, for all a G N'^, 1 < |a| < A^, we have, 


We assume also that /x satisfies (18). 





e 

2fi 


(l + eV|VC|W 


= eV ■ (VV^(„) + e5“CVw) - -wd^G^{xf) 

h 

-/3Xwd%!f^^{dtb) + Sc,. 


Furthermore Sa is controlled 


r OO 


\^S^\l 2 < eMNS^iU) + C (^Mn, ^ e£^iU)"2 + Mn 15^61^0 

Proof. The proof of this Proposition is similar to the proof of Proposition 4.10 in m 


expect we use Propositions B.8 and B.13 See also Proposition 6.4 in [16] . 


□ 


Thanks to this linearization, we can ” quasilinarize” equations (13). It is the purpose of 


the next proposition. Let us introduce, the Rayleigh-Taylor coefficient 


a := a{U, /3b) =1 -h edt ( w[eC, fib] { ip, —dtb 


+ £^V[eC, /3b] ( V’, ^9tb\ • V TwfeC, /3b] 


This quantity plays an important role. We also introduce two new operators, 

A[U,/3b] : 


0 -j^G,[eC,pb] 
a{U,/3b) 0 


and 


B[U,pb] := 


eV • (.P) 0 

0 eP-V 


(25) 


(26) 


(27) 


We can now quasilinearize the water waves equations. We use the same arguments as in 
Proposition 4.10 in |I9] and part 6 in [T6] . Notice that we give here a precise estimate 
with respect to dtb and P of the residuals Ra and Sa and that the minimal value of N, 
regularity of U, is smaller than in Proposition 4.10 in |19j . 


Proposition 2.8. Let T > 0, N > max (to, 1) + 3, 6 G P G 

L°°(]R+; i7'^+^(]R‘^)) and U G satisfies (19) for all 0 < t < T and solving (13). IPe 


assume also that /i satisfies (18). Then, for all a G N“, 1 < |a| < N, we have. 
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t 


dtU^^)+A[U,f3b]iUia)) + l{ic.i=N}mmUia)) = [ ^Gf^^[0,0Wdtb),-d<^P 


+ ( Ra, Sc 


(28) 


Furthermore, Ra and Sa satisfy 


Ra\L^ < MN\{eC,l3b)\HN£^{Uy2 + —Mm \dtb\L^M^ > 


/?A 


|^ 5 a|i 2 < eMMS^iU) + C^N,Y + Mm ■ 


Proof. Thanks to Proposition B.15, we get 


\G^^[eC,/3b]{d°‘dtb) - G^^[0,0]{d‘^dtb)\^^< / \dGf^^[z£C, zl3b]{d°‘dtb).{C,b)\^2dz 

J 0 

< MM\{£C,,ldb)\HN \d°‘dtb\^^HN . 

Then, denoting Ra = Ra + G^^[eC,, I3b]{d°‘dtb) — Gj^^[0,0]{d°‘dtb), we obtain the first 


equation thanks to Proposition |2.6[ For the second equation, using Proposition 2.7 and 
the first equation of the water waves problem, we have 


= -a“C - ev ■ (yf;^a) + ed^CVw) + +Pwd'^G^^{dtb) -d^P +Sa 

h- 

= -(9“C - £V ■ (VV’(a) + ed^fVw) + swdtd^C - d^P + Sa 

= —(9"C(1 + sdtw + • Vw) — eV • V'0(a) + £dt{'wd°‘C,) — d°‘P + Sa 

= -a5“C - sV ■ VV^(„) + edtiwd^C) - d^P + 5«, 

and the result follows. 

□ 

In the case of a constant pressure at the surface and a fixed bottom, it is well-known 
that system (28) is symmetrizable if 

3 Clmm ^ 0 , a{U,/3b) > 

Clmin* (29) 

Then, we introduce the symmetrizer 
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This symmetrization has an associated energy 


jr-{U) = - (5[C/,/36](C/(„)),C/(„))^,, if a ^ 0, 

= kiS + ^ Uh, ^G^[eC,Pb]iAh)) , 

i H'i i \ fj. y ^2 

j^W(U)= 

\a\<N 


(31) 


As in Lemma 4.27 in [19], it can be shown that and are equivalent in the 
following sense. 


Proposition 2.9. LetT > 0, N £ N, U £ Eif satisfying (19) and (29) for all 0 <t<T. 
Then, for all 0 < k < N, is comparable to 


1 


< £^{U) < (mn + —) F^'^^[U,b]. (32) 

\ ^min J 


\a{U,mL^ + Mn 

2.4 Local existence 

The water water equations can be written as follow : 

dtU + N{u) = (0,-p)*, 

with M{U) = iAfi{U),Af 2 {U)Y and 

M,{U) := -^G,[eC,/3bm - ^G^^[eC, mdtb), 


(33) 


Af 2 {U) := C + ^ (l + e^^|VCH (^w[eC,/3b] 


(34) 


According to our quasilinearization, we need that a be a positive real number. Therefore, 
we have to express a without partial derivative with respect to t, particularly when t = 0. 
It is easy to check that (we adopt the notation of Remark B.12 in Appendix |B.3|) 


aiU, /3b) = 1+ e"R[eC, /3b] ( f;, ^dtb ] • V 


w[£C,Pb] { V’, ^dtb 


+ £dw[i/j,^dtl^ . {-Afi{U),dtb) + ew[£C,/3b] (-P - Af2{U), ^dfb 


(35) 


The following Proposition gives estimates for a{U,/3b). It is adapted from Proposition 
6.6 in IHI. 


Proposition 2.10. Let T>0, to >^,N> max(to, 1) + 3, (C, V’) ^ ® solution of 

the water waves equations (13), P £ L°°(]R+; i7'^+^(]R‘^)) and b £ 1T^’°°(M+; i7'^( 
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such that Condition (19) is satisfied. We assume also that // satisfies (18). Then, for 
all 0 <t <T, 


\a{U,/3b) — llffto <c(MN,m.ax{/3X, (3) \dtb\j^^jjN,eS^{U)' 2 j £S^{U )‘2 
+ £Mn(\VP\ ■ 

Furthermore, if dfb G L“(M+; and dtP G L°°{R+; {R<^)), then, 


\dtia{U,f3b))\Hto <C[MM,ma^iPX,fi)\dtb\^i,^j,^,\VP\L^HM,££^\U)-^]££^\U)-2 

\ t X t X / 

+ eC'(MAr,max(/3A,/3) \dtb\^^jjN)(\S/P\,^i,^j^N+— \d'fb\^i,ooj^X 


Proof. Using the first point of Proposition B.8 and Product estimate C.2 we have 


\V_[£C„fib]{£fi,f5Xdtb)X[w[£C,,l3b\ {£'if,l3Xdth)]\^to<MN +fiX\dtb\^^^t^ 


Furthermore, thanks to the hrst point of Proposition B.15 and the first point of Theorem 
3.28 in [19] we obtain 


edw{fi,—dtb].{-Afi{U),dtb) 


<Mn\ {eNi {U),fidtb) I jyto+i \dtb\ 




Then, the first inequality follows easily from Proposition |B.8[ Proposition |2.5| and Prod¬ 
uct estimate C.2 The second inequality can be proved similarly. □ 


We can now prove Theorems 2.3 and 2.4 We recall that S := max(e, /3^). 


Proof. We slice up this proof in three parts. First we regularize and symmetrize the 
equations, then we find some energy estimates and finally we conclude by convergence. 
We only give the energy estimates in this paper and a carefully study of the nonlinearities 
of the water waves equations is done. We refer to the proof of Theorem 4.16 in [19] 
for the regularization, the convergence and the uniqueness (see also part 7 in |16)). 


For Theorem 2.3 




respectively on 
[0, T] ( respectively on 


respectively Theorem 2.4), we assume that U solves (13) on [0, T] 


^ and that (19) and (29) are satisfied for ^nd on 
for some T > 0. 




a) |a| = 0, The 0 - energy 


We proceed as in Subsection 4.3.4.3 in m and part 6 in [16]. We have 
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- (Al{M2iU)-C),-G,[eC,l3b]i^h)) - (-G,[eC,f3b]i^h),^"^p) • 

V d- /L2 7^2 

(36) 

We have to control all the term in the r.h.s. 


* Control of ^ (^A^G^^leC, (3b]{dtb), A 2 (^'j 
Using Proposition |B.7| we get 


L2 




AlG^^[eC,midtb),Ak)^^ 


< MN—\dtb\^ooHN £^{U)l 

S t X 


* Control of (A 2 (A/ 2 (t/) - C), hGt^l^C, 


L2 


Using Proposition 2.5 and Proposition B.8 we get 


AHM2{U)-C),^G^[eC,Pb]{Ah) 

d- /l2 


<\M2iU)-C\^3 


^Gf,[eC,Pb]{A2^l;) 

A 


L2 


<eMjvf^(U)i+ Mn \dtb\L^Hj^^ e£^{U). 


* Control of ^^G^[eC,/36](A2?/)), AaP^^^. 
We get, using Remark 3.13 in m, 


-G^K,/3^](a4),A^^’ 

/L2 


< MnE^{U)2\VP\l^jjN. 


* Control of ^ (^dGfj,[eC, pb]{A 2 'ijj).{dtC,dtb), A 2 'ij;^^^. 


Using Proposition 3.29 in |19j . the second point of Theorem 3.15 in [19], Proposition B.7 


and Proposition 2.5 , we get 


^ (dG^ [eC,f3b]{A2'ijj).{dtC,dtb),A2 


< MM\{eMi{U),/3dtb)\HN-2 


<Mfqe£ {U )2 +max(/3,/3A) \dtb\i^^^N £ (U). 
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Finally, gathering all the previous estimates, we get that 


< sMNS^iuf^ +MnC (p^ax, \dtb\L^H^) I3)S^{U) 


(37) 


b) |a| >0, the higher orders energies 

We proceed as in Subsection 4.3.4.3 in [19] and part 6 in m- A simple computation 
gives 


= -el{|a|=w} (aCH, V • (ChZ))^. + ( aC(„), ^G^^[0,0]{dtd‘^b) + R. 


L2 


-el{|a|=iv} f-G;,[eC,/36](V’(a)), [Z-VV'(a)]) + (-Gf,[eC, Pb]{'ip(^^)), Sa - d^P 

VM /l2 \fJ- 7^2 


(38) 


We have to control all the term in the r.h.s. 
* Control of (9toC(„), C(a))i 2 - 


Using the second point of Proposition 2.10 we get 


(9iaC(„),C(a))^2 < M^G Lax, , \VP\^^j,^,sS^{U)"A ef^([/)i 


+ C{MN,l3X\dtb\j^^jjN^i^\S/P\yyi,oc^N+^ \d^b\^i,o<,^N'j eS^{U). 


■k Control of (^aC(a), ^Gf^^[0,0]{dtd°b)^ 


We get, thanks to Proposition 2.10 and B.7 


^ <C'(5max,/Umax,|&| 


1\/3A, 


{Up. 


-k Controls of el{|„|= 7 v} (aC(a), V • (C(a)Z)L- 

Inspired by Subsection 4.3.4.3 in [I^, a simple computation gives 
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^(^C(ck) (C(Q;)i^) )^2 ^ (^C(a)^ * C{a)^j^2 

<c(^Prne..,f^n...,\b\^y2,^H^,\VP\L^HN,SS^iU)'j£[S^{U)l+S^{U) 
using Proposition |2.10] and Proposition |B 


* Controls of (^G/,[eC,/3^](V’(a)),5’a-5“p)^^, (^^dG/,[eC,/36](V>(a)).(«9tC,V’(a))^ 2 ^nd 


L2 


Rr 


L2 


We can use the same arguments as in the third and the fourth point of part a) using 
Propositions |2.8| and |2.10[ 

* Control of e (^^G^[eC,/3b](V'(a)), [Z • VV’(a)])^ 2 - 

We refer to the Subsection 4.3.4.3 and Proposition 3.30 in m for this control. 


Gathering the previous estimates and using Proposition 2.9 we obtain that 


_d 

dt 


^ 'linin CXmin t X t X / 


•*min 

eF^{U) 5 +max(e, P)F\U)+P^{U) $ 




(39) 

Then, we easily prove Theorem |2.3t using the same arguments as Subsection 4.3.4.4 in 
[19]. Furthermore, for a G [O, ^], defining {U){t) = 5‘^°‘R^{U) (^), we get 


^ R^{U) < C'fpn,ax,l^max, — ,^, 161^3,00^^ , ,R^iU) 

' Omin n ■ X y''t 


dr 


mm "'mm 


We can also apply the same arguments as Subsection 4.3.4.4 in m and Theorem |2.4| 
follows. □ 


2.5 Hamiltonian system 


In this section we prove that the water waves problem (13) is a Hamiltonian system 


in the Sobolev framework. This extends the classical result of Zakharov (|33|) to the 
case where the bottom is moving and the atmospheric pressure is not constant (see also 
|in|l. In the case of a moving bottom, P. Guyenne and D. P. Nicholls already pointed 
out it in naD We have to introduce the Dirichlet-Dirichlet and the Neumann-Dirichlet 
operators 


^It seems that there is a typo in their hamiltonian; ’’—(v” should read 
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( 40 ) 


\ GneC./3'.l(W = ('I''*)|„_.+^t, 

where is defined in © and is defined in p^ . We postpone the study of these 
operators to appendix [B| 

Remark 2.11. If we denote := <f>'^ + ^ $ satisfies 


2 ^ = 0 in ; 


= V’ , y/l + W\Wdn%=-i+^b = —dtb. 


Then 


and 


^l+e‘^\VC\^dn%=e^ = + —G^^^[ecmdtb)^ (41) 


= G^^H^fibm + ^G^;^[eCmdtb). 


r^ND r 


(42) 


Theorem 2.12. LetT > 0, to > 5 , C, ^ e C°([0, T]; R*o+i(M'^)), fj e 

dtb G C°([0, r]; R^(M'^)), P G C'^([0, T]; L^(]R'^)) such that (CjV’) is a solution of (13). 

Define H = H{(,Tp) = T(C,V’) +^(C;^); where T{C,'fi) = T is 


T= 


1 


‘^hJQt 
and U{Q^ fi) =U is 


7^4 / fU'S' _L 


V£J ch- + 


+ 




dtb [G^^[eC,fibm + ^G^;!^[eC,mdtb)j, 

(43) 

(44) 


11 = - CdX + / (PdX. 

^ .md ./rod 


Then, the water waves equations (13) take the form 


dt 


cUfo 

\-i 0) [d^Hj ■ 

Remark 2.13. T is the sum of the kinetic energy and the moving bottom contribution 
and lA the sum of the potential energy and the pressure contribution. Using Green’s 


formula and Remark 2.11 we obtain that 


T=- 


1 f ,/l„r. -- fix r 


-G^[eC,fibm + l^GrK,/36](9t6) dX 


+ 11 —dtb (Gj^^[eC,fibm + ^G^^[eC,fib]idtb) ) dX, 

2 . Kd e \ ^ e ^ 
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Proof. Using the linearity of the Dirichlet-Neumann and the Dirichlet-Dirichlet opera¬ 
tors with respect to if and the fact that the adjoint of G^^[eC, fib] is fib] (see 

Proposition |B.5|), we get that 


d^H = ^G^[eC,l3bm + ^G^^[eC,mdtb). 


Applying Proposition B.13 (which provides explicit expressions for shape derivatives) 
and remark 2.13, we obtain that 


2d(^H = ^b]{f^)w + • U - e^G^^[eC, pb]{dtb)w + 2P + 2C, 

= -[eC, /36] {i^)w + eVV’ • VV’ - [sf, pb] {dtb)w + 2P + 2C, 

= e IVV^I^ - -w^ (1 + eVlVCH +2P + 2C, 

which ends the proof. □ 

In fact, working in the Beppo Levi framework for ^f requires that T^\^dtb G L^(M‘^) and 
results that are not dealing with this paper. 


3 Asymptotic models 

In this part, we derive some asymptotic models in order to model two different types 
of tsunamis. The most important phenomenon that we want to catch is the Proudman 
resonance (see for instance j24| or m for an explanation of the Proudman resonance) and 
the submarine landslide tsunami phenomenon (see [20] . |28] or [29]). These resonances 
occur in a linear case. The duration of the resonance depends on the phenomenon. For a 
meteotsunami, the duration of the resonance corresponds to the time the meteorological 
disturbance takes to reach the coast (see |24jL However, for a landslide tsunami, the 
duration of the resonance corresponds to the duration of the landslide (which depends 
on the size of the slope, see m or [28]). If the landslide is offshore, it is unreasonable 
to assume that the duration of the landslide is the time the water waves take to reach 
the coast. A variation of the pressure of 1 hPa creates a water wave of 1 cm whereas a 
moving bottom of 1 cm tends to create a water wave of 1 cm. Therefore we assume in 
the following that abott,m = a (and hence fiX = e). However, it is important to notice 
that even if for storms, a variation of the pressure of 100 hPa is very huge, it is quite 
ordinary that a submarine landslide have a thickness of 1 m. Typically, a storm makes 
a variation of few Hpa, and the thickness of a submarine landslide is few dm (we refer 
to |20)). In this part, we only study the propagation of such phenomena. Therefore, 
we take d = 1. In the following, we give three linear asymptotic models of the water 
waves equations and we give examples of pressures and moving bottoms that create a 
resonance. The pressure at the surface P and the moving bottom bm move from the 
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left to the right. We consider that the system is initially at rest. We start this part by 
giving an asymptotic expansion with respect to /r and max(e,/3) of Ph]. 


Proposition 3.1. Letto > i, and b G such that Condition (19) is satisfied. 


We suppose that the parameters e, fi and pL satisfy (18). Then, for all B ^ H 
with 0 < s < to + I; we have 




\G^^[eC,f5b]{B)-G^^m{B)\^^_. <Mo\{eC, 1^1 

and 


Proof. The first inequality follows from Proposition |B.15| and the second from Remark 

[O 

□ 

Remark 3.2. In the same way and under the assumptions of the previous proposition, 
we can prove that (see Proposition 3.28 in m), for 0 < s < to + I; 

|GmK,/3^](V’)-G^[O,O]('0)|^,_i < ^Mo|(eC,/36)|iyio+2|Wl^s+i 

and 


-Gi,[0,0](fi) + Afi 
h 


^<piG\Vf;\ 


5 . 

2 


We denote by V the vertically averaged horizontal component, 


V = V[eC,fib]i^|;,^tb) 


1 

1 + e(^ — fib 



Vx{<^[£C,fib] {fi,dtb) {■,z))dz, 


where = ^[sC,fib] {if,dtb) satisfies 


(45) 


j - ^ + fib< z <£(, 

\ = V' , yi + fi‘^\Vbfidn%=-i+/3b = hdtb. 

The following Proposition is Remark 3.36 and a small adaptation of Proposition 3.37 
and Lemma 5.4 in m (see also Subsection A.5.5 in m)- 


Proposition 3.3. Let T > 0, to > 5 , 0 < s < to and C,b G ([ 0 , T]; R*o+2(M'^)) 

such that Condition (19) is satisfied on [0,T]. ITe suppose that the parameters £, fi and 
/i satisfy ([Ts]). We also assume that fit G W^'°° ("[0, T]; Lf®+^(M'^)y Then, 


G^[£C,fibm + piG^^^[£C,fib]{dth) = -pTJ-{{l + eC-fib)V)+pidtb, 
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and 


1^ - - ^C'f)^>/^max,e|CliT‘0+2,/3|fe|^oo^^0+2j max (^\V'iIj\^s +2 , \dtb\];^^jj^+i^ 


min 

dtV— IJ,C )/^maX)|CI_H'*0+2 j^to+2 ,\b\y^^2,oo^^+2 , |V'0| j^s+2 , \dt^'>P\}{s+2 

In this part, we will consider symmetrizable linear hyperbolic systems of the first order. 
We refer to [7] for more details about the wellposedness. In the following, we will only 
give the energy associated to the symmetrization. 


3.1 A shallow water model when (3 is small 

3.1.1 Linear asymptotic 

We consider the case that e, (5, ^ are small. Physically, this means that we consider 
small amplitudes for the surface and the bottom (compared to the mean depth) and 
waves with large wavelengths (compared to the mean depth). The asymptotic regime 
(in the sense of Definition 4.19 in [19]) is 


Aw = {(e, /3, A, /x), 0 < ia,£, P < 6o, /3A = e}, 


(46) 


with ho <C 1. 


Proposition 3.4. Let to>i, N> max(l,to) + 3, t/° G E^, P gW^’^{R+;H^+\R<^)) 
and b G W^’°°{R'^; {R‘^)). We suppose (19) and (29) are satisfied initially. Then, 

there exists T > 0, such that for all {e,f3,X,fj.) G Alw, there exists a solution U = 
(C, Ip) G to the water waves equations with initial data and this solution is 

unique. Furthermore, for all a G [O, |), 


c-c 

/ 


x + 

Vip — Vip 

/ 


\ 


L°° ( 

n PI1_ 

D, SOL 

L 0 J 

;jyiV-4(]Rd) j 


L°°( 

n — 

r^o J 

.J^iV-2(Rd) J 




where 


1 1 


C — C Is (t/°), , , |6|^3 .oo^^, |VP|^i,oo^^ j , 

\ ^min t a t a y 

and with, (C, ip) solution of the waves equation 


with initial data . 


dtC + Axip = dtb, 
dtfi + C = -P, 


(47) 
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Proof. First, the system (47) is wellposed since it can be symmetrized thanks to the 
energy 


£{t) 



+ 




2 

L2 ■ 


2.4 


Using Theorem 
equation and for a. 


we get a uniform time of existence > 0 for the water waves 
1 parameters in Alw- Then, using Proposition 3.1 Remark 3.2 


Proposition B.8 and C.l and standard controls we get that 


with 


dtC + AxV’ = dtb + Ri, 

+ C = -P + P 2 , 


(48) 


j \Ri\hn-4 < c(^e\C\HN,\b\^^HNy\{eC, ^b)\HN + ^) max , 

\ \R2\hn--l < eC(e|C|H^, |6|i^^v)max 

If we denote Ci = C “ C and tpi = tp — ip, we see that (Ci) V’l) satisfies 


Differentiating the energy 


dtCi + ‘^x'tpi = Ri, 

dtipi + Ci = R2- 


£^{t) = l\Ci\lN-, + l\Vip,\l^., 


we get the estimate thanks to Proposition |2.5| and energy estimate in Theorem 2.4 □ 

This model is well-known in the physics literature (see j25|l. 


3.1.2 Resonance in shallow waters when f3 is small 


We consider the equation (47) for d 
equation for h := (p — b, 


1. We transform it in order to have a unique 


' d'^h - d\h = d\{P + b), 

< /i|t=o =-6(0,.), (49) 

, dth\t=Q = 0. 

We denote f{t,X) := {P + b) {t,X), which represents a disturbance. We want to un¬ 
derstand the resonance for landslide and meteo tsunamis. In both cases, it is a linear 
respond, in the shallow water case, of a body of water due to a moving pressure or a 
moving bottom, when the speed of the storm or the landslide is close to the typical wave 
celerity (here 1). We can compute h thanks to the d’Alembert’s formula 
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h{t, X) = -^ (6(0, X -t) + 6(0, X + t)) + ]^ f dxfir, X + t- T)dT 

V J 0 


hT{t,X) 


■=hL{t,X) 


\ [ dxf{T,X -t + T)dT. 

^ Jo 


■=hR{t,X) 

We are interesting in disturbances / moving from the left to the right (propagation to a 
coast). Therefore, we study only hn- The following Proposition shows that a disturbance 
moving with a speed equal to 1 makes appear a resonance. 

Proposition 3.5. Let f G L°°(M+; 

X eR, t > 0, 


and dxf G ^ R‘^). Then, for all 


\hR{t,X)\<^\dxfL. 

Furthermore, if f{t,X) = fo{X -t), fo G H^{R^) and |/o(^o - io)| = |/'loo equality 
holds for {to,XQ). If f{t,X) = fo{X - Ut) with fo G and [7/1, 

<min(^i:^^,^|/'|^^ . 

Proof. lff{t,X) = fo{X-Ut), 

hnit, X) = fo{X - t + (1 - U)T)dT, 

and the result follows. □ 

This Proposition corresponds to the historical work of J. Proudman ([25]). We rediscover 
the fact that the resonance occurs if the speed of the disturbance is 1. For a disturbance 
with a speed different from 1, we notice a saturation effect (also pointed out in [28]). 
The graph in Figuregives the typical evolution of \h{t,-)\^ with respect to the time 
t for different values of the speed. We can see the saturation effect. We compute h 
with a finite difference method and we take f{t,X) = . We see also that 

the landslide resonance and the Proudman resonance have the same effects. There are 
however two important differences that we exposed in the introduction of this part. The 
first one is the duration of the resonance. A landslide is quicker than a meteorological 
effect. The second one, is the fact that the typical size of the landslide (few dm) is bigger 
than the size of a storm (few hPa). For instance, for a moving storm which creates a 
variation of the pressure of 3 hPa during 15toj the final wave can reach a amplitude of 
13 cm (it is for example the case of the meteotsunami in Nagasaki in 1979, see |24j). 
Conversely, an offshore landslide with a thickness of 1 m that lasts to, can create a wave 
of 50 cm (which corresponds to the results in [28|1. Therefore, we see that the principal 
difference between an offshore landslide and a moving storm is the size. 
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Figure 1: Evolution of the maximum of h, solution of equation (49), with different values 
of the speed U. 


3.2 A shallow water model when /3 is large 
3.2.1 Linear asymptotic 

In this case, we suppose only that s and ^ are small. We recall that /36(t, X) = j3bo{X) + 
PXbmit, X). Then, we assume also that I — bo > > 0. In the following, we denote 

ho := 1 — l3bo. The asymptotic regime is 


Alvw = {(e, /3, \, fi), 0 < e, < 6o, 0 < P < 1, (3\ = e} , 
with do 1. We can now give a asymptotic model. 


(50) 


Proposition 3.6. Let to > ^, N > max(l,to) + 4, 6 G = 

(Co)'0o) £ Eq, and P We suppose that (19) and (29) are sat¬ 

isfied initially. We suppose also that bo G {W^) and that /iq = 1 — jdbo > hmin- Then, 
there exists T > 0, such that for all (e,/I,A,/i) G Alvw, there exists a unique solution 
U = (C) V') £ water waves equations with initial data U^. Furthermore, for V 

as in (45), 


1C CllL°o([0,r];H^-4(]Rd)) + < ^<50(5, 


where 


C = c(£" (UO) I VF| 

\ /^min ^min t a t a y 

and (Cl, El) solution of the waves equation 
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( 51 ) 


' dtCl + V • [h^Vi) = dtbm 
< + VCi =-VP, 


_ (Ci)|t=o = Co, (F)|f=o = V [sCo,Pb\t=o] (V’o, (9i&)|t=o) 


Proof. The system (51) is wellposed since it can be symmetrized thanks to the energy 


For the inequality, we proceed as in Proposition |3.4[ differentiating the energy 


£^{t) = I \C2\lN-, + ^ (hoA^-^V2,A^-^V2)^, , 


with ((2 = C — Cl and V 2 = V — Vi- 
standard controls, we get result. 


Using Gronwall’s Lemma, Proposition 


3.3 


and 

□ 


This model is well-known in the physics literature to investigate the landslide tsunami 
phenomenon (see [28] 1. 


3.2.2 Amplification in shallow waters when j3 is large 

In this part, d = 1 and we suppose that P = 0. The same study can be done for a non 
constant pressure. For the sake of simplicity, we assume also that initially the velocity 
of the landslide is zero and hence that (9t5m)|t=o ~ ^ bottom does not move at the 
beginning). We transform the system © in order to get an equation for only. We 
obtain that satisfies 


dtCi - dx (do^xCi) = dfbm, (52) 

with (Ci)|t=o = 0 and (<9tCi)|j=o ~ wonder now if we can catch an elevation of the 

sea level with this asymptotic model. Therefore, we are looking for solutions of the form 

C2{t,X) = tC3{t,X). (53) 

The following proposition gives example of such solutions for bounded moving bottoms 
(with finite energy). 

Proposition 3.7. Suppose that ho > hmin > 0 with ho e P^(M). Let be a 

solution of 

( dtCz + dx {hoV3) = 0, 

\ dtV 3 + dxC,‘s = 0 , 

with ^ P^(M). Then, Ci{t,X) = tC, 3 {t,X) is a non trivial 

solution of (|52|) with 
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bm{t,X) = 2 C 3 {s,X)ds, 


( 54 ) 


and bm{t, •) is bounded in and in L“(M'^) uniformly with respect to t 


\bm{ti ■)!L°° — 


where C is independent on t. 


Proof. Plugging the expression of f and bm in (52), we get the first result. We have to 
show that Cs £ Consider the linear hyperbolic equation 


I dtrj + dx {h^W) = 0, 

[ dtW + dxv = 0 , 

with (r?, hP)|i=o ~ (“/)0). This system has a unique solution (r],W) G (R] (R)). 

Furthermore, {dt'r],dtW) G L^(M)), and {dtrj^dtW) satisfies the same linear hyper¬ 

bolic system as (Csj^s)- By uniqueness, Cs = dtd and 


Since, for all t. 


b^{t,X) = 2p{t,X) + 2f{X). 


r]{t, Xf + ho{X)W{t, XfdX 


fiXfdX, 


and ho > /imin > 0, we get the control of OIlz- Finally, p satisfies the waves 

equation 


dh - dx (hodxv) = 0 , 

with (f/,9 j?7)|^^q = (—/,0) G H^{R'^). Then, for all t, 


\dtT]{t,Xf + ho{X) \dxv{t,X)fdX 


ho{X)f{XfdX. 


Therefore, \i]\jji (and |r/|^oo by Sobolev embedding) is controlled uniformly with respect 
to t. □ 


In the following, we compute numerically some solutions of Equations ( |52[ ) of the form 
( [^ with a finite difference method. We take bo{X) = — tanh(X), /? = | and (i9tC3)|t=o ~ 
(4X^ — 2)e“^^. The figure I is the evolution of the maximum of • The figure is 
the graph at different times of the waves and the landslide. The dashed curves are the 
landslide, the solid curves are the waves and the dotted curve is the slope. Therefore, 
we see that an important elevation of the sea level is possible even if we do not consider 
that the seabed is flat. 
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Remark 3.8. In order to simplify, we consider that the system is initially at rest. But 
our .study can easily be extended to waves with non trivial initial data. In particular, we 
can study a wave amplified by a landslide. This is what happened during the tsunami in 
Fukushima in 2011 (see We compute numerically this amplification. We consider a 

wave moving with a speed equal to 1 (typical speed in the sea after nondimensionalization) 
that is amplified by a landslide. Figure represents the evolution of the maximum of 
this wave. We can see an amplification. 



Figure 2: Evolution of the maximum of Ci, solution of (52), for a non flat bottom bo. 



Figure 3: Evolution of the surface Ci (solid line), solution of (52), and the landslide 
(dashed line). 



Figure 4: Evolution of the maximum of h, solution of (52), with non trivial initial data 
and with bm like in Figure]^ 
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3.3 Linear asymptotic and resonance in intermediate depths 

In this case, we consider only that e, /3 are small. Physically, this means that we consider 
small amplitudes for the surface and the bottom (compared to the mean depth) and that 
the depth is not small compared to wavelength of the waves. In this part, we generalize 
the Proudman resonance in deeper waters. The asymptotic regime is 


■^LWW = {(e,/3, A,/r), 0 < e, /3 < (5o, /3A = e and 0 < ^ < ^max} , (55) 

with Ao <C 1 and 0 < /imax- Using the energy 


S{t) = \ \C\h + \ , 

and proceeding as in Proposition |3.4| (we need also Proposition 3.12 in |19j). we get a 
new asymptotic model. 


Proposition 3.9. Let to > ^, N > max(l,fo) + 3, 6 G VU^’“(M+; (M*^)), [7° = 

{Co,i’o) £ and P G IU^’“(M+;iif'^+^(M'^)). IPe suppose that (19) and (29) are 
satisfied initially. Then, there exists T > 0, such that for all {e, j3, p) ^ Alww, there 
exists a unique solution U = (C, fi’) G E\ to the water waves equations with initial data 

U^. Furthermore, for all a G [O, 




1^1 ( / 

c-c 




'*0 


.HN-2 




where 


C = cU^ (17°) , —, ^max, |VP| 

V ^min l^min * ^ 


W, 


jjN 


where 


is a solution of the waves equation 


1 


dtC--G'M[0,0](V^) = G;:^[0,0](di6), 
M 

dtfi + C = -P, 


with initial data . 


(56) 


The Proudman resonance is a phenomenon which occurs in shallow water regime. We 
wonder if there is also a resonance in deeper waters. In this part, we only work with a 
non constant pressure and hence dfi) = 0. The same study can be done for a moving 
bottom. We consider the equation (56) for d = 1. Since, the initial data does not affect 
the possible resonance, we suppose in the following that = 0. We transform the 


29 

















system (56) in order to have a unique equation for (in the following we denote C by C 
to simplify the notation) 

—G’mIOjOKC) = —G^[0,0](P), 

H /r 

C|t=o = 0) dtC\t=o = 0- 

We can solve explicitly the previous equation, we get that 


««) = ,/f 






dr 


:=a(i,0 






dr. 


:=Ci?(iiO 


In order to find a resonant pressure, we suppose that P has the form e where a is 

a real smooth odd function which is sublinear, there exists C > 0 such that |a(^)| < Cl^l. 
We also suppose that the phase velocity of the disturbance is positive, > 0. Pq is 

a smooth function in a Sobolev space with Po{0) / 0. We denote ^ 

simple computation gives that 


Furthermore, we have 


ia(t,-)i < ia(t,-)iLi < 


Pn 


P 


lCR(i,OI = 2 

t 

< 

“ 2 


with an equality if and only if a(^) = ^a;(y^^). Hence, it is natural to consider that 


P{t,^) = (57) 

A simple computation gives 

CRit,X) = -| / (58) 

^ Jr 

We wonder now if a resonance occurs. We need a dispersion estimate for the linear water 
waves equation. 
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Proposition 3.10. Let f £ such that /(O) = 0. Then, 





< 


Vt 



1 


+ //8 


Li(R) 


1^1^ 



Proof. We denote I{t), 


I{t) := f 

JR 


1 

y/h 


[ e-^My‘^^y^-Ty)f(JL)dy. 

Jr \y/Tj 


We denote (f>, 


(t^iv) = v^{y) - —V, 


and ?/o the unique minimum of (j)". Figure [^represents (j)" on [0, +oo[. 



Figure 5: Profile of 4>" ■ 


To estimate I{t) we decompose I{t) into four parts. 


h{t) 


VT^Jo \VtJ 


-/ 

\/h Jo 

I rvo 

h Jo 


d 

dy 


rvo 


e ^ 'dz 


' ivt) 
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Then, using Van der Corput’s Lemma (see m) and the fact that for z & [y, yo] 
\(t)"{z)\ > W'{y)\ and \ct)"{z)\ > Cz, 


c j_ /-y / 


lAWI < . r ! 

Jo 

c /■+“ 

< 


dy 


y/V-Vi Jo 

Furthermore, for M > y^ large enough. 


y/^if) (6 


dc 


rM 


hit]=J=r 

y/y- Jyo 


e f ( ^ ] dy 


fM 


r-M 


A( / e 

dy 


y_ 

y/i^J 


'yo 


1 /-A^ ry 

1 / ! —1 


Jyo y/J^ V y/T^ J Jyo Jyo 


e 


4>{z) 


dz 



dy. 


Then, using Van der Corput’s Lemma and the fact that for z £ [yo,y], 

\r{z)\>\cP"iy)\stnd\riz)\>Cz-l, 



Tending M to +oo we get the result. The control for .^ < 0 is similar. □ 

Therefore, in the linear case, we have also a resonance. 

Corollary 3.11. Let Pq £ n such that XPq £ H^(K) and let 

0 < y < ^max- Consider, 

Cnit,X) = ~ [ 
z Jr 

Then, 

ICfl(i) Oloo ^ C'(/^max)y^ {\Po\h3 + I-^oIlI + l-^-Poli^s) , 

and 


lim 

t->+oo 




> C{Po) > 0. 
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Proof. We take f{f) = Then, 



<{l + ^\f\) PoiO +kl 



and the first inequality follows from the previous Proposition. For the second inequality, 
we use a stationary phase approximation. We denote (/>(^) = Let > 0, such 

that ^oTb(Co) = ^Pq ^ , and < 0, such that Then, we have. 


L°° 


lim 

t->+CXD 


1 

Vt 


CR{t,tX^f) 


lim ^ 

t^+co 2fl 


2/i4 


^ujiOPo 
^{f,0^/Jj-)CoPo{^o) 


A_ 




VWWM 


Since |<?i''(OI ^ T'|^| and a;(^ov^) > C'(?o)\//^> result. 


□ 


Remark 3.12. Notice that for all s G M, 

— |VTo|j|^s + 2 . 

Hence, by tending formally [i to 0, we rediscover the result we get in the shallow water 
case (section \3. l[ j. 

Remark 3.13. Notice that for a general pressure term P{t,X) we can show that the 
amplitude Q satisfying 


Cnit, •) + 2^o(' ~ 




tanh(^|g|) ^^^^ 




dr 




^^K\ 


dr, 


satisfies also 


K(L OIoo — (|T’lLoo(R+.^i(Rd)) + |T’|2,oo(]K+;H3(R‘*) 


+ 




Hence, contrary to the shallow water case, we can not hope a linear amplification with 
respect to the time t. Corollary 3.11 also shows that the factor of amplification of \fi is 
optimal. 
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Hence, we observe that in intermediate water depths, a resonance can occur but with a 
factor of amplification of \/t and not t. But we saw that in the shallow water case, the 
resonance occurs for a moving pressure with a speed equal to 1, P{t,X) = Po{X — t). 
We wonder if this pressure can create a resonance. The following Proposition shows that 
the previous pressure can create a resonance with a factor of amplification of fs. 

Proposition 3.14. Let 0 < < //max- Let Pq G L^(M) n such that Po{0) / 0. 

Consider, the amplitude (r created by P{t,X) = Po{X — t), 


Ut ,0 = I 


l-t 


Then, 


Oloo — C'(/tmax) f — l-Polii + I-PoIhI^ • 


Furthermore, if XPq G P[^{R), 


(59) 


Proof. We have 


lim 

—Cnit, •) 

C 

> — 

MO) 

t-^+oo 

t3 

CO 



^ ./tr> 


I-t 




= _/ ^ujP^Pq —^ e I e"^ 

h Jr V \/T‘- J J-t 


e ^ dsd^. 


We decompose this integral into 3 parts. 


\hit)\ = 


1 


1 

t3 

< — 




^ fioiOPo ( e [ e" ^^e'^^d^ds 


Vt 


i-t 


Furthermore, since |w(^) — 1| > for 0 < |^| < 1, we have 


|/2(t)| = 


hJt-h<\^\<i \y/hj J-t 


< C 


1. [ Pq ( 

1 

^3 


Pn 
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Finally, 


l^3(t)| = 


J\i\>i Vv^y J-t 


1 

^ J\i\> 

1 ^ ^(0 


e 


y/J^J\^\>l w( 0 -l V\/^ 


Po ( ^ ) fe - e ) d^ 


< C 




PoiO dC, 


— C'//4 [Pol/fl ) 

and the first inequality follows. For the second inequality, we use a stationary phase 
approximation. Wedenote0(^) := ^(a;(^) —1). We recall that (^(^) = —g.^^+o(.^^). Using 
a generalization of Morse Lemma at the order 3, there exists a > 0 and V' £ C°° ([—a, a]), 
such that for all |y| < a. 


Then, 


4>{'^{y)) = ^0"'(O)y^,'(/'(O) = 0 and V’'(0) = 1. 
b 




f ^'(2/)w(V'(y))V’(y)Po 

J —a 






e dy + o{s a) 


= (hll\ " Po(o) [ zP^'dz + o(s-i). 
\ S J 


Therefore, 


^o(O) 


Then, in intermediate water depths, a traveling pressure with a constant speed equal to 1 
is also resonant, but it takes more time to obtain a significant elevation of the level of the 
sea. In the following, we compute numerically some solutions. We take Po(y^) = —e~^ 
and fj. = 1. The figure is the evolution of a water wave because of a pressure of the 
form ( [ST] ). The solid curve is the wave and the dashed curve is the moving pressure. 
The figure is the evolution is a water wave when the pressure moves with a speed 1. 
The figure compares the evolution of the maximum of the resonant case and the case 
when the speed is equal to 1. 


lim 

t->+oo 


\CR{t,t) 

t3 
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Figure 6: Evolution of the surface elevation in 
moving pressure P in (57) (dashed line). 


(solid line) because of a resonant 



Figure 7: Evolution of the surface elevation in (59) (solid line) because of a moving 
pressure P with a speed of 1 (dashed line). 



Figure 8: Evolution of the maximum of Cj? in the resonant case (solid line) and the 
moving pressure with a speed of 1 (dashed line). 


Remark 3.15. In our work, we neglect the Coriolis effect. However, in view of the 
duration of the meteotsunami phenomenon, it would be more realistic to consider it. It 
will he studied in a future work (J2^) based on the work of A. Castro and D. Cannes 
and 
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A The Laplace problem 

A.l Formulation of the problems 

In this part, we recall some results of Chapter 2 in m and Section 4 of [16] and study 
the Laplace problem (© in the Beppo Levi spaces. We suppose that the parameters e, 
IX and /3 satisfy Condition ( [I^ . The Laplace problem ( [ItI ) is 

f Af = 0 in SI, , 

\ = 0, yiTswan-i-®i„-i+s, = B, 

where B = dth. Notice that is here the upward conormal derivative 




'^X,z^^\z=-1+Pb- 


For the study of (15) we refer to m- We work with Beppo Levi spaces. We refer to 
m and Proposition 2.3 in |19j for general results about these spaces. We recall that, 
for s > 0, 


:= E VV' E , 

that H^{R'^ X (—1,0))/]^ is a Hilbert space for the norm \Vx,z • \l^ and that H^{R'^) is 
dense in H^{R'^). In order to fix the domain, we transform the problem into variable 
coefficients elliptic problem on S := R^ x (—1,0) (the flat strip). We introduce a reg¬ 
ularizing diffeomorphism. Let 0 : M ^ M be a positive, compactly supported, smooth, 
even function equal to one near 0. For <5 > 0 we define 

5 _ 

{X,z)^{X,z + a{X,z)), 

and 


cj(A, z) := [e{Sz\D\)eC{X) - 9{6{z + l)\D\)/3b{X)] z + e9{dz\D\)C{X). 

We omit the dependence on t here. In the following, we denote by M a constant of the 
form 


7 ) /^max, £|Clj7*0+l(Rrf)) /5|^li7*0+l(E‘i) 

^min 

In order to study the Laplace problems in S, we have to treat the regularity in the 
direction X and in the direction z one at a time. We introduce the following spaces. 

Definition A.l. Let s E M. We define (H*’^(S'), Ms,!) and Ho,!) 

H^’\S) := LlH^xiS) C and 
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and 


:= LlH^xiS), and \u\j^s,o = 

Remark A.2. We have the following embedding (see Proposition 2.10 in mJforsGR 

H^+h\s) C LfH^iS). 

In the following, we fix 5 > 0 small enough. Then, we can transform our equations. We 
denote := and we get that 


with P{T‘) 


r = 0 in S, 

\ 4>%=0 = 0 , dn4>^\z=-l — 
Id+ixd+1 + <3(S) and 


(60) 


Q(S):= 


Idxd 

-yJJiVxcr^ 


-y/Jl^xcr 

-dzcr+^i\'Vx<^\^ 

l+dzcr 


(61) 


Notice that T’(S) is well defined if 6 is small enough and that cln := e^ • (P(E)V^^ • ). 
We have to know the regularity of T’(S). It is the subject of the next proposition (see 
Proposition 2.18 and Lemma 2.26 in |19jL 


Proposition A.3. Let to > 

Then, 


i, C^b £ i/*“+^(]R‘^) such that Condition (19) is satisfied. 




(s) 


AS) 


Furthermore, P(S) is coercive. There exist a constant A:(S) > 0 such that < M and 

V0 E , V(A, z)£S , P(S)(A, z)0 • 0 > fc(S)|0|2. 

We have a variational formulation of the Laplace problem ( |60[ ). We introduce 

R1.../(5) := V{SC{z = -l}) ' = V{SC{z = -l}) ' . 

See Proposition 2.3 (3) in [19] for a proof of the second equality. 


Definition A.4. Let B £ H 2 ( 


We say that 4> £ H^surfi^) ® variational 


solution of (60) if for all (p £ Hq surf(^)’ 


We have also the following trace result that we can prove easily using a density argument. 
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Lemma A.5. For all (p G we have 


\Jl + v^l^l F\z=-i 


< 2 


L 2 (R‘^) 




L^S) 


We can now establish existence and uniqueness results. 


Proposition A.6. Let B ^ H 2 (M*^) and C, b £ satisfying 

problem (60) has a unique variational solution named surfiS). 


Then, the 


Proof. Because S is bounded in the direction z and that P(S) is uniformly coercive, the 
results follow from the Lax-Milgram theorem and Poincare inequality in Hq □ 


In this part, we study the Laplace problem ( 17|) 
(see Chapter 2 in m) and we can transform 


, but the same work can be done for (15) 
(17) as follows 


= 0 in 5, 

(l)%=o = = 0 . 


In the following, we denote by , the unique solution of (62). 


(62) 


A.2 Regularity estimates of the solutions 

In this part, we give some regularity estimates. 

Theorem A.7. Let to > ^ and 0 < s < to + 5 - Let C,b £ be such that 

Condition (19) is satisfied. Then, for all B £ we have 

1 




L2(S) 


< M 


V^l + ^/Ji\D\ 


B 




Futhermore, if s > max(0,1 — to); 'we have 



< M 

_ 1 _ B 


LHS) 

V^ + ^\D\ 




Proof. Let d > 0 and y be a smooth compactly supported real function that is equal 
to 1 near 0. We introduce the smoothing operator A| := x((iA)A®. We know that 
B^ £ Rq^^^j(S). Therefore, using A^^B^ a test function, we have 


[ V^^B^ • PiJ^M^^AfB^ = - [ B{AfB%^_,. 

Js URd- 

Since P(S) is symmetric, A| commutes with ^ and is independent of z we obtain 
that 
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A 1 


jRd -^1 + y^|-D| 

Then by coercivity of P(S) and trace inequality |A.5 


bUi + ^\d\aib^ 


b=-i 


fc(s) 




L 2 (S) 


< 


+ 2 


[A|,g(s)]v£,i?‘ 

d X,2 ^2(5) 


1,2 






yr+V^ 


L2(5) 

B , 

L 2 


and 


MS) 




< 


L^S) 

We have to distinguish two cases. 

a) 0 < s < to : 


[AlQ{^)]Vl^B^ 


L^S) 


+ 2 


Af 


xA+VM 


B 


L 2 


The commutator estimate C .6 (with Tq = to) and Proposition A.3 give 


MS) 


^sK,.B^ 


L2(5) 


<^^IQA)L 


=M°(s) 




< M 




LHS) 


+ 2 


L^S) 

Ai 


+ 2 




v^i+VAA 


B 


L 2 


\/i + \AAI 


B 


L 2 


for some e > 0 small enough (e < to — 5 ). Using a finite induction on s and taking the 
limit when 6 goes to 0, the first inequality follows. For the second estimate, we only 
need to give a control of d^B^. We use Equation (60) satisfied by B^. We express P(S) 
as 



{1 + a{X,z))Idxd 

q\X,z) 


A simple computation gives 


q(X, z) \ 
1 + Qd+i{X, z)J 


(1 + qd+i)d^,B^ = - ^Vx • ((1 + a)^XxB^) - ^Vx • {d,B\) 

- ^/JldzCi ■ XxB^ - y/Jid^XxB^ • q - d^qd+idzB°. 
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We have a, q, qd+i G LfH^^{S), d^qd+i G LfH^^ ^{S) and 1 + g^+i > K'^)- 

Then, since s > 1 — to and Vx-B® G by the product estimates C.3 and C.4 (with 

Tq = Iq), we obtain the result. 


b) to < s < to + 5 


The commutator estimate 


C.7 


(with To = to + 5 and ti > ^) and Proposition 


A.3 


give 


MS) 




LHS) 


< M 




LHS) 


+ 


A 




+ 2 


LHS) 


A| 


y^T+7W\ 


B 


L2 


We denote e ■.= — ti. We obtain the first inequality for to < s < to + e thanks to the 

previous case. Furthermore, we saw that 


(1 + qd+i)dlB^ = - • ((1 + a)^VxB^) - ^Vx • (9.B®q) 

- \//^<9;jq • VxB® - • q - d^qd+id^B^. 

We have a, q, qd+i G {S), ct^q, dzqd+i G ^ (S) and 1 + qd+i > k 

s-i r 

Then, since s > 1 — to and VxB® G ^ (S), by the product estimates C.3 and 

(with To = to), and we obtain the second inequality for to < s < to + e. Using a finite 
induction, we obtain the first and the second inequality. □ 


s). 


C.5 


B The Dirichlet-Neumann and the Neumann-Neumann op¬ 
erators 

We refers to Chapter 3 of m for more details about the Dirichlet-Neumann operator 
and Section 3 in m for the study of these operators. 


B.l Main properties 

We can express the Neumann-Neumann operator with the formalism of the previous 
section. For tp ^ (M'^) and B £ (M'^) we have 

G^K,/36](^) = (e.-P(S)Vj^,,V’'’), ^ 

\ ’ / | 2:=0 


(63) 


and 


Gfi^[eC,Pb]iB) = (^ez ■ PimlzB"') 


\z=0 


(64) 
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Remark B.l. Notice that (see Proposition 3.2 in m) 


iG,[0,01W = ^ CriO, 01(B) = 


B. 


We recall that G^[eC,^ /35] is symmetric and maps continuously H 2 (M'^) into yH 2 (F'^'' 

(see Paragraph 3.1. in USD- We need an extension result in H 2 (F'^) in order to give 
a dual formulation of the Neumann-Neumann operator. 


Definition B.2. Let ip £ H' 


We define as 


# _ sinh([z + 1]^|D|) 
sinh(^|D|) 


Remark B.3. satisfies weakly 

I = 0znS, 

I P%=o = P*\z=-i = 0 - 

We can prove easily regularity results for similar to . 
Proposition B.4. Let s > 0 and (p G Then, 




1 

H- — 

LfiS) y/Jl 




# 


LfiS) 


< C 


I + ./Jl\D\ip 


Hs-i 


We can now give a dual formulation of the Neumann-Neumann operator. We introduce 
the Dirichlet-Dirichlet operator, for fit £ H 2 


— Lifi 


h=-i 


Gl!^[eC,fibm ■= 

The following result is Proposition 3.3 in [El- 

Proposition B.5. Let to > ^, B £ H~ 

satisfied. G^^[£(, fib]{-) can be extended to H~ 2 {W^) with the dual formulation 


(65) 


and C,b £ such that (19) is 


G^^[eC,mB)={ 


H2\ 






( 66 ) 


Furthermore, the adjoint of G^^ [eC, fib] is Gjf^[eC,fib]. For all B £ H 2 (F“) and 
ip £ Ri(F‘^), 

{GfH,fib]{B),ip) . . = {B,G^^H,fib]{ip)) , . . 
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In order to study shape derivatives of the Dirichlet-Neumann and the Neumann-Neumann 
operators, we have to introduce the Neumann-Dirichlet operator. For B G we 

define 


(67) 

The following result is a symmetry property and a dual formulation of the Neumann- 
Dirichlet operator. 


Proposition B. 6 . Let B £ H 

G^^[eC, pb]{B) can be view as 


and C-,b £ 77*°+^ (M'^) such that (19) is satisfied. 


G^^[eC,mB) 




( 68 ) 


Furthermore, G^^ [eC, fib ](•) is a negative symmetric operator and, for all Bi, B 2 in 

i7-i(M‘^), 


{G^^^[eC,fib]{B^),B2)^^_,,,y_^_,,, = (G;y^K,/3&](i?2),i?i)(^_v2y_^-i/2 • 

We refer to Proposition 3.3 in M for a proof of this result. 

B.2 Regularity Estimates 

In this part we give some controls the Neumann-Neumann operators. 

Proposition B.7. Let to > 5 ; 0 < s < to + 5 o.'nd C,,b £ such that Condition 


(19) is satisfied. Then, Gj^^[£C,fib] maps continuously 2 

\G^^[eC,fibm\^^_^ <M\B 


into itself 




Proof. This Proposition follows by Theorem A.7 and by using the same arguments that 
Theorem 3.15 in [19]. □ 

We can extend these estimates to w[£C,fib], the vertical velocity at the surface and 
to P[eC) fib] the horizontal velocity at the surface. These operators appear naturally 
when we differentiate the Dirichlet-Neumann and the Neumann-Neumann operator with 
respect to the surface C- We dehne 


w[£C,fib] ■■= 






>-£■ 


{R^)^ 77"-2 (M'^) 

G^[eC, fibm + fib]{B) + EfiVC • (69) 

l + eV|VCP 
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and 




(V’,5) ^V'ip-sw[eC,f3b]{^l;,B)VC. 


(70) 


Proposition B.8. Let to > i, 0 < s < to + I and £ H^°~^^(R'^) such that Con¬ 


dition (19) is satisfied. Then, w[£C,l3b] maps continuously 2 \ 


H^-kR‘^' 


X H^-2\ 


into 


\w[eC,,fib]{fi,B)\^^_i < M + lA 


B\ 




Furthermore, */1 < s < to; maps continuously H^~^^{R'^) x 2 (M*^) into 


H^-2 


1 + IBI 1 


< M/x 

Finally, we have the same continuity result for V_[eC,, fib]. 

We can also give some regularity estimates for G^^[£C,, fib] since it is the adjoint of 

G^,^H:fih]. 

Proposition B.9. Let to > 5; 0 < s < to + 5 o,nd C,b € such that Condition 

(19) is satisfied. Then, GjA^[£(,fib] maps eontinuously H^~^ 2 {R'^) into itself 


\VG^^[eC,fibm\^^_r < M|VxA|^,_i . 

Finally, we can give some regularity estimates for G^^leC, fib]. 

Proposition B.IO. Let to > |, 0 < s < to + | andC,,b G such that Condition 

(19) is satisfied. Then, GjA^[e(,fib] maps continuously H^~ 2 (R'^) into 


In the same way, we can extend also these estimates to wfeC,, fib], the vertical velocity at 
the bottom and to P[eC, fib] the horizontal velocity at the bottom. These operators ap¬ 
pear naturally when we differentiate the Dirichlet-Neumann and the Neumann-Neumann 
operator with respect to the bottom b 

TB + fipSIb-V{G^^H,fihm + yiG^^[eC,fib]{B)) 
w[£C„fib]{'fi,B) = - , , n-) lvTi .19 -’ (^1) 


1 + ,5V|V6|2 


and 


YXeC.fibfifi’^B) = V {G]A^[£C,fibm + ^iG;^H,fib][B)) - fiw[£C,fib]{fi^,B)Vb. (72) 
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Proposition B.ll. Let to > f; 0 < s < to + 5 CL'^d C^b € such that Con¬ 
dition (19) is satisfied. Then, ud^C^fib] maps continuously x into 


\w[eC,fib]{'ip,B)\^^_i < M • 

Finally, we have the same continuity result for VleC, fib]. 


B.3 Shape derivatives 

Let to > Given B & H 2 We denote by L the set of functions (C, b) in 


satisfying (19). We introduce the map 


G^iB) := 


H2{ 

-iNNt 


(73) 


iC,b)e^G^;^^[eC,fib]{B), 

which is the Neumann-Neumann operator. We can also define G^{'ip), w{'il),B) and 

V{d,B). 

Remark B.12. When no confusion is possible and to the sake of simplicity, we write 
Gjd^iB), w{i^,B) andV{'fi,B) instead of G ^[eC, fibfiif), Gjd^[eC,fib]{B), 
w[eC, fib^ip, B) and VfeC, fififiip, B). 

In order to linearize the water waves equations, we need a shape derivative formula for the 
Dirichlet-Neumann and the Neumann-Neumann operators. The following proposition is 
a summarize of Theorems 3.5 and 3.6 in 1161 and Theorem 3.21 in 


Then, 


Proposition B.13. Let to > 'l,C,b€ € Hi{W^) and B e 

G^(V’) and G^^{B) are Frechet differentiable. For {h,k) G we have 


dG^ifi).{h, 0) + fidG^d 0) = - fib]{hw[eC, fibfifi, B)) 

-sfiV ■{hV[eC,fib]{if,B)), 


and 


dG^ifi).{0,k)+frdGf^^{B).{0,k) = fiiaGfd^[sC,fib] (v ■ [k V[eC, fibfifi, B))) . 
Furthermore, 


dGf{if).{h,t)) + p.dG^^^{B).{h,t)) = -eG^^H,fib]{hw[eC,fih]{if,B)). 

Thanks to these formulae we can give some controls to the hrst shape derivatives of the 
operators. For instance, we give an estimate for dpi and dY_. 
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Proposition B.14. Let to > 5 o,nd {C,b) G such that Condition (19) is 

satisfied. Then, for 0 < s < to + for fi G (M'^) and B G we have 


dV{fi,B).{h,k) 1 , \dw{'tlj,B).{h,k)\ i< M\{h,k)\fjto+i 

H 2 -W 




Proof. This result follows from Proposition B.13 and Proposition B.7 


□ 


We end this part by giving some controls of the shape derivatives of and . We 
do not use the previous method, we differentiate j times directly the dual formulation 
of both operators. We refer to Proposition 3.28 in [19] for a control of d^G^.ih., k)(^/!). 


Proposition B.15. Let to > 2 (Cj^) ^ such that Condition (19) is 

satisfied. Then for all 0 < s < to + ^ and B G (1^^)? we have 


\d^G^^.{h,k){B)\^^_r<Ml[ I (.^hi , \ jjtQ+i 

i>l 


\B\ 




Furthermore, if 0 < s < to and B G 77*° 


i>2 


We do not prove this Proposition here (which is based on a shape derivative of B®). We 
refer to 


C Useful Estimates 


In this part, we give some useful estimates, product and commutator estimates. We 

refer to Appendix B in |19] . |18j and Chapter II in for the proofs. The first estimates 

are useful to control ^/. We recall that fC = 

^ Vi+vTl^l 


Proposition C.l. Let f G 


and g € H2 


Then, 


\^g\L 2 <g IWI^i <max(l,/i 4 )|V/|i ,2 and |V/|i 2 < max(l, ^ 4 )|rp/|^i. 

Proof. The first inequality follows from the fact that 1 + the second 

inequality from < max(l, and the third from ^ max(l,/x 4 ). □ 

We need some product estimates in W^. The following Proposition is Proposition 2.1.2 

in |5]. 
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Proposition C.2. Let s,si,S 2 £ ^ such that s < si, s < S 2 , si + S 2 > 0 and s < 
Si + S 2 — §■ Then, there exists a constant C > 0 such that for all f G and for 

all g G we have fg G and 


Ifgln^ < C\f\H‘>i\g\H‘’2- 

We also need some product estimates in S := x (—1,0). The following Proposition 
is the Corollary B.5 in m- 

Proposition C.3. Let s, si, S 2 G M such that s < si, s < S 2 , si + S 2 > s < si + S 2 — ^ 
and p G {2, +oo}. Then, there exists a constant C > 0 such that for all f G L'^Hff{S) 
and for all g G L^H^^{S), we have fg G L^H^{S) and 


1^'^ {fa) Il?L2_(S) ^ 

The following propositions gives estimates for \/{l-\- g) in the flat strip S. We refer to 

Corollary B.6 in 1191- 

Proposition C.4. Let Tq > —To < s < Tq, ko > 0 and p G {2, +oo}. Then, for all 
f G LfzHx{S) and g G L^Hjf{S) with I + g > ko, we have 


A* 


/ 


1+5 


LW^{S) 






Proposition C.5. Let Tq > 2 , s > —Tq and ko > 0. 
H^’^{S) with 1 + g > ko, we have 


^ Wf\LU^{sy 

Then, for all f,g G LfH^°{S) n 


/ 

1+5 



+ l-{s>ro}l/l 


LP-Hi 


Tolaln^’O 


Notice that if s < To, f G H^’^{S) is enough. 


We need some commutator estimates in S. The following Propositions are Corollary 
B.17 in HH]. 

Proposition C.6. LetTo > <5 > 0, 0 < ti < 1 withti < Tq-| and< s< To + ti. 

Then for all u G and v G we have 


|[A|, 


u 


v\l2(s) ^ ^ 


Afu 


L^Lj^iS) 


|A: 


s—ti 


Proposition C.7. Let To > |, (5 > 0, 0 < ti < 1 with ti < To 
Then for all u G ffAj.o y g we have 


L^S) ■ 

— I and 


^ < S < Tq + ti. 


|[A|,tt] v\j^2(^s) < C 


Afu 


L^S) 


I A' 


S — ti 


^LfiS) 
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